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We review the latest variational calculations of the ground state properties of doubly closed shell 
^\ • nuclei, from ^^C to ^"^Pb, with semirealistic and realistic two- and three-nucleon interactions. The 

O^ I studies are carried on within the framework of the correlated basis function theory and integral 

^^ ■ equations technique, with state dependent correlations having central and tensor components. We 

report results for the ground state energy, one- and two-body densities and static structure func- 
tions. For ^®0 and ''"Ca we use modern interactions and find that the accuracy of the method is 
O . comparable to that attained in nuclear matter with similar hamiltonians, giving nuclei underbound 

\^ ' by ~2 MeV/A. The computed Coulomb sums are in complete agreement with the latest analysis of 

^ , the experimental data. 
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I. INTRODUCTION: INTERACTION AND CORRELATIONS 



Our knowledge of the nuclear interaction is steadily improving because of both the ever increasing number of 
experimental data and the more and more sophisticated theories which are being developed to tackle the longstanding 
^5 problem of an accurate description of the strongly interacting nuclear systems. The interplay between experiments 
y—{ ' and theory is clear: interactions are built which fit the data (typically from nucleon-nucleon scattering) and are 
^\ , then tested theoretically in more complicated structures. This approach has led to the construction of nucleon- 
0^ ' nucleon (NN) potentials which reproduce a huge amount of NN scattering data (the latest versions fit '^1800 pp and 
~2500 np data with x^ ~1 ^'^d break the charge independence and charge symmetry ||l|-^). However, the use of 
these accurate, modern NN potentials in A>2 nuclei has also proved their inability to correctly describe the nuclear 
binding. Light nuclei are underbound and the nuclear matter saturation density is not correctly reproduced by a 
hamiltonian containing two-nucleon forces only. Several cures to these pathologies have been proposed, ranging from 
relativistic effects to extra degrees of freedom (A's) or to many-body forces. 

If we choose to advocate the lack of binding to the presence of more-nucleon potentials, then the first step consists 
in the introduction of three-nucleon interactions (TNI). The knowledge of TNI is by far less deep than that of the NN 
potential because of the considerably lower number of experimental data and of the higher difficulty in constructing 
^ ■ theoretical models. However, even with these caveats, the present models of TNI provide the correct binding for 
Cu ^ A=3,4 nuclei (actually, TNI contain parameters which are fitted on light nuclei binding energy) and, more important 
and less obvious, bring the computed nuclear matter saturation density very close to its empirical value of pnm =0.16 
fm-3. 

Following these guidelines, one can write a realistic nuclear hamiltonian in the form: 



i<j i<j<k 



2m 
where Vij and Vijk are two- and three-nucleon potentials. The NN potential may be split into: 



V^J=V1^+V:,+V^ (2) 

where the electromagnetic part, ff" contains Coulomb, Darwin-Foldy, vacuum polarization and magnetic moments 
contributions; the long range part is well established from a microscopic point of view as being given by the One Pion 
Exchange potential (OPEP), vf^; the intermediate and short range potential, v(j, arises from the exchange of more 
pions and heavier mesons (p and uj), however its microscopic derivation from meson exchange theory is troublesome 
(for the larger number of Feynman diagrams involved and the greater difficulty in the potential extraction) and 
a phenomenological approach is often followed, consisting in a physically plausible parametrization in terms of a 



reasonable number of parameters. For instance, the Argonne vis NN potential Q| contains 43 parameters and it is 
given by a sum of 18 terms 

^ii = Yl ^pi^^3)0fj , (3) 

p=l,18 

where the first 14 terms are isoscalar, with 

Of=i'i4 ^ [1^ a^ . a,,S,j,L ■ S, L2, lV,; • a,, (L • S)^] ® [1, r, ■ r,] , (4) 

and the remanining isovector and isotensor components are: 

Of~^^'^^ = [l,<Ji-aj,Sij][3Ti^Tj^ -Ti-Tj],[Ti:,+Tj^] . (5) 

Two-pion exchange gives large part of the attraction of the three-nucleon potential, whereas the remainder is 
mostly phenomenological. Recently, attempts are being made to evaluate microscopically other pieces, as the three- 
pion two-Delta diagram Q . A popular form of TNI are the Urbana forces |5| , given by 



«*jfe=<fe+<fe ■ (6) 



The two-pion Fujita-Miyazawa Q term, w?'^., is 



"§*: = X! ^27r{Xij,Xik}{Ti ■ Tj , T^ ■ Tfe } + C27r [Xij,Xik] [t^ • Tj,Ti ■ Tfc] , (7) 

cycl 



with 



Xij = Y^{rij)a, ■ a J + T^{r.^)Sij , (8) 



y,r(r) and T^(r) being the Yukawa and tensor Yukawa functions. The short range part, vf,^^ is written as 



v}Jk=Y.UoTl{n,)n{nk) ■ (9) 

cycl 

^2ir, C2-n 6 C/q are chosen to reproduce the binding of A=3,4 nuclei. 

Once a realistic hamiltonian has been built, task of the theorist is to address the ground state Schroedinger equation, 
H'^Q — i^o^Oj to solve it or to find good approximations to its solutions. The exact solution is now possible in light 
nuclei by a variety of techniques: Green's Function Monte Carlo 0| for 3<A<8, Faddeev and Faddeev-Yakubovsky 
g and Correlated Hyperspherical expansions (CHE) 0. GFMC has been pushed to the highest A-value, whereas the 
other approaches have been used mainly for A=3,4. However, they are more flexible in the sense that both Faddeev 
and Correlated Hyperspherical expansion methods are more easibly extendible to the study of three- and four- 
body scattering reactions. In fact, several few-body reactions of relevant astrophysical interest have been accurately 
analyzed by CHE. For the remaining part of the nuclear table, including nuclear matter, methods to exactly solve 
the Schroedinger equation are still to come, even if a very promising approach, based on the Hubbard-Stratonovich 
transformation, is currently under development [|lO| . Hence, alternative ways to gather information on the many-body 
wave function must be devised. The Brueckner theory has brought standard perturbative theories to an extremely 
high level of sophistication ||ll|]; the Coupled Cluster (CC) method has been used to build correlations into the wave 
function and is, by now, a standard technique in many-body physics [Oj; the variational principle has provided a 
powerful recipe to construct accurate wave functions. 

A direct consequence of the strong nuclear interaction is the failure of the mean field (MF) approach in the 
description of the nucleus. A striking evidence of this fact is provided by (e, e')p experiments |T3] where a clear 
signature of a depletion of the occupation probability with respect to unity for states below the Fermi level, ep, (as 
requested by MF theories) is present, together with a non zero occupation probability for states above ep- This 
behavior may be explained in terms of short range dynamical correlations (in contrast with statistical correlations, 
due to the antisymmetry) generated in the wave functions by the interaction and which can be hardly described by 
standard perturbation methods based on a non interacting basis. For instance, Brueckner theory must sum infinite 
numbers of ladder diagrams in order to deal efhciently with the nuclear potential. 

In the Correlated Basis Function (CBF) theory the non perturbative correlation effects are directly embedded into 
the basis functions. This property makes the theory a powerful tool to investigate many-body interacting systems in 
several fields of physics, as liquid Helium, electronic structures (both in the forms of electron fluids and lattices) and 



both finite nuclei and infinite nuclear matter. The flexibility of the CBF approach results in a realistic description 
not only of the ground state (energy, momentum distribution, distribution functions and so on) but also of dynamical 
(cross sections) quantities. 

A set of correlated basis wave functions, ^n(l, 2. ..A), may be built by applying a many-body correlation operator, 
F{1,2...A), to the model basis functions, <i>„(l, 2...^), 

*„(1,2...A)=F(1,2...^)$„(1,2...A) , (10) 

where the operator F is intended to take care of the dynamical correlations, whereas the model wave functions, <&„, 
include antisymmetrization effects and, possibly, long range correlations due to collective excitationsfas BCS type 
states or surface effects). A perturbative theory may ve developed in terms of the correlated states (p^), having the 
nice property of a rapid convergence and the counterindication of a large difficulty in computing the matrix elements. 
The zeroth order of this theory is often referred to as the variational level of CBF. In fact, the variational principle 
is used to fix the correlation operator by minimizing the ground state energy, Eq = (^o|-ff|^o)/(*o|*o)- 

The choice of the correlation operator depends, to a large extent, on the interaction. A form of F{1, 2. ..A) that has 
shown to be suitable to nuclear systems is 

F{l,2...A)=s\ n P^n ' (11) 

i.e. a symmetrized product of two-body correlation operators, F^. The model states, $„(1,2...A), are Slater de- 
terminants of single particle wave functions, (f>aii), obtained by some MF potential (simple plane waves in infinite, 
homogeneous matter) . Fij is chosen, consistently with the interaction, to be given by 

F^, - E f'(^v)0!, , (12) 

p=l,8 

where the sum runs up to the spin-orbit components in eq.(0). If the p > 2 components are disregarded, the Jastrow 
scalar correlation is recovered. The free minimization of Eq would provide the best choice for the correlation functions 
f^(r). However, this is a prohibitive task in nuclear systems, so the correlation functions are usually parametrized 
and the parameters are fixed by minimizing Eq . An additional minimization may be performed on the parameters of 
the MF model wave function. 

II. SOMETHING ON CLUSTER EXPANSIONS AND FERMI HYPERNETTED CHAIN EQUATIONS 

As already mentioned in the previous Section, computing the matrix elements in the CBF approach represents the 
greatest difficulty in its application. In principle one could use Monte Carlo (MC) based algorithms to sample the 
occuring multi-dimensional integrals. For simple interactions and correlations this can be done even for very large 
numbers of particles. In such a case, the A = oo limit can also be studied by means of appropriate conditions at the 
borders of the simulation box. Liquid Helium, both as a fluid and as droplets, has been subject of MC investigations. 
However, the strong state dependence of the nuclear interaction prevents from the use of MC techniques in medium- 
heavy nuclei, as well as in nuclear matter. The ^^O ground state was studied within CBF by the cluster Monte 
Carlo (CMC) method [|lj], where the Jastrow correlations contribution is exactly treated by MC sampling, and that 
from the remaining operatorial components is approximated by considering (via MC) up to four- or five-body cluster 
terms. 

An alternative approach is provided by cluster expansions. The expectation values of operators are written in terms 
of n-body densities, 

p^il)^{J2S{r.-rl)) , (13) 



p^i\l,2) = {J26{v,-r,)Sir,-r2)0!^) , (14) 

and related quantities, as the density matrices. The densities are then cluster expanded in terms of dynamical, 
h{r) = [/^(r-)]^ — 1 and f^{r)fP~^{r), and statistical, po{i,j) = J2a ^a(*)'?^a(i)' correlations. The expansion is linked, 



in the sense that disconnected diagrams, given by the product of two or more pieces not connected among each other 
by any kind of correlation and coming from the expansion of the numerator, cancel exactly with those generated by 
the expansion of the denominator. However, the cluster expansion is not irreducible, i.e. the surviving connected 
diagrams bear vertex corrections (except for the Jastrow correlation in the A=oo case, where the vertex correction is 
simply the density). 

The cluster diagrams are then classified into: Nodal, N, where all the paths between the external points (for 
instance, points 1 and 2 for the two-body density of eq.(l4)) go through the same internal point, or node; Composite, 
C, given by the product of two or more nodals; the remaining diagrams belong to the Elementary class, E. 

Fermi Hypernetted Chain (FHNC) integral equations [ISl allow for the exact summation of the nodal and composite 
diagrams, in the Jastrow case, once the sum of the elementary diagrams, E{1, 2), is given. So, the function E{1, 2) is 
actually an input for solving thr FHNC equations. Unfortunately, no exact way of computing this function is known 
and one must use approximations. The simplest one is the FHNC/0 truncation, where the choice E = is made. 
This seemingly crude approximation is however appropriate to nuclear systems, where the density is not very high. 
In denser liquid Helium, a more accurate treatment of the elementary diagrams is needed. 

For state dependent correlations, a complete FHNC summation is possible only for the Jastrow part. In fact, 
different orderings of the operators in the symmetrized product (O) may give the same cluster contribution and a 
scheme to correctly keep track of all of them has not been devised so far. However, partial classes of diagrams containing 
operatorial correlations may be exactly summed by the Single Operator Chain approximation (FHNC/SOC) |1^. The 
FHNC/SOC integral equations sum all the nodal diagrams with only one operatorial correlation per internal side, 
besides all the Jastrow correlated clusters. The accuracy of the FHNC/SOC approximation has been set to ^1 MeV/A 
in nuclear matter at saturation density by computing its leading corrections ||l^. The FHNC/SOC scheme has been 
recently extended to doubly closed shell nuclei in Is coupling p^ . 

A way to control the approximations is provided by the density sum rules. In fact, the exact densities satisfy the 
following normalizations: 

fdripiil)^A , (15) 

^ ^dn I dv2p^^\l,2) ^ 1 , (16) 



A{A-l) 



^jdv,jdvj^\l,2) = -l . (17) 

The last equality holds for isospin saturated systems. Deviations from the sum rules are a powerful check for the 
approximations made in the solution of the FHNC equations, both in the treatment of the elementary diagrams and 
of the operatorial correlations. 

As a matter of fact, FHNC/SOC satisfies the normalizations of the densities to a high degree of accuracy, as it 
is shown in Table |. The Table gives the one-body (Si, eq. (|l5|))and central and isospin two-body (S2 and Sr, 
eqs. ( p^Jl7| ) density normalizations for Jastrow (/i) and operatorial (/g) correlations in ^^O and ^°Ca. The /g model 
does not contain spin-orbit components but it does have the tensor ones. The Jastrow model satisfies practically 
exactly the sum rules, whereas the SOC approximation produces discrepancies of ~9% in the worst case |lq| . 

III. ENERGIES, DENSITIES AND STRUCTURE FUNCTIONS 



In a series of papers [119|-pi[| we extended the FHNC formalism to doubly closed shell nuclei, both in Is and jj 
coupling scheme, described by Jastrow correlated wave functions. In those papers model interactions were used, as 



the semirealistics Afnan and Tang (S3) one |22|. The original S3 potential is purely central, without tensor components, 
and reproduces the NN s-waves scattering data up to ~60 MeV. Its modified version |2^ has been supplemented in 
the odd channels of the same repulsion as in the even ones and gives a nuclear matter binding curve close to the current 
realistic potentials. The correlations we adopted contained a r^ dependence, allowing for distinguishing between the 
different nucleon pairs (/„„ 7^ fpp ^ fnp)- The ground state energies for the Tz-dependent correlations are shown in 
Table O together with the experimental values. The correlations have been fixed by solving the Euler equations which 
result from the minimization of the expectation value of the hamiltonian at the second order of the cluster expansion. 
Given the poor quality of the interaction, the comparison with the experiment is not really meaningful, even if not 
completely unsatisfactory. However, it is interesting to note a saturation trend of the binding energy along the mass 



number. The Table also gives the energies obtained by an isospin independent correlation, identical for all the pairs 
(Average Correlation Approximation, AC A). The poorer variational quality of the ACA correlation produces higher 
energies than the isospin dependent model. The density normalizations are always satisfied with an accuracy greater 
than percent for the calculations of the Table. 

000000000 

FHNC/SOC has been pushed to deal with realistic potentials and state dependent, non central correlations in ^^O 
and '^'^Ca in Refs. ^8 24|. In the last work a complete hamiltonian of the form given in eq.(||) with the inclusion of 



the momentum dependent parts of the NN interaction and of the three-nucleon potential, as well as of the Coulomb 



term, was adopted. Table III gives the g.s. energies per nucleon for two models: Argonne vu + Urbana VII [ES 26 
(A14+UVII) and Argonne Wg + Urbana IX 0] (A8'+UIX). The correlated wave function features a /e nuclear matter 
correlation and single particle wave functions generated by a Woods-Saxon mean field potential. In the A14+UVII 
case, the ^^O FHNC/SOC energies are compared with the CMC results of Ref. ||l^, obtained with a similar wave 
function, and with the Coupled Cluster ones of Ref. |p^. The r.m.s. radius is also given. The variational A14+UVII 
energies are close to the CC estimates, while the discrepancies with CMC are due to the approximations in the 
FHNC/SOC scheme and to small differencies in the correlation. Actually, our nuclear matter Euler equation docs 
not contain the spin-orbit potential, whereas that used for the CMC correlations does. This fact produces small 
variations, which mostly affect the kinetic energy. FHNC and CMC give similar radii, both of them smaller than CC. 
Explicit three-body correlations, not given by the product of two-body ones, were found to provide ^^0.85 MeV/A 
extra binding @ in the A14+UVII case. 

The charge densities, Pc{ri), and the two-body distribution functions, P2(?'i2), computed with the A8'+UIX wave 
function are shown in Figure^. The CBF Pc(''i) (dashed lines) are folded with the proton form factor and compared 
with the experimental data (solid lines). /92(^i2) is defined as 

P2{ri2) = jld'Ri2pi'\r,,r2) , (18) 

where R12 is the center of mass coordinate. Short range correlations strongly affect P2{ri2) (solid lines) at small 
internucleon distances, where the NN repulsion heavily depletes the distribution functions with respect to the IPM 
estimates (dashed lines). 

Two-body densities are also of great interest since they allow for analyzing several integrated nuclear cross sections. 
In fact, the responses of a nucleus to an external probe, described by an operator Ox , can be expressed in terms of the 
dynamical structure functions, Sx{q,^), whose non energy weighted sum gives the static structure function (SSF), 
Sxiq)- In turn, the SSF are given by expectation values of appropriate combinations of one- and two-body densities. 
For instance, the electromagnetic longitudinal response, 5L(q,aj), as measured in electron-nucleus experiments, is 
mostly due to charge fluctuations. The corresponding SSF, ^^(q) (or Coulomb sum), is given by 



SUq) = 1 + - / d^n I d^r2 e''!--- [ppp(ri, r2) - Pciri)Pcir2)] , (19) 



where Ppp (ri,r2) is the proton-proton two-body density. 

In Figure g we compare the A8'+UIX CBF Coulomb sums with those extracted from the world data on inclusive 
quasi-elastic electron scattering |28] experiments in -'^^C, ^"Ca, and ^^Fe. The Figure also shows the nuclear matter 
Coulomb sum for the A14+UVII model from Ref. [g9|. The agreement of the nuclei theoretical SSF with the latest 
experimental data, which properly take into account the large energy tail of the response, is complete, whereas the 
nuclear matter shows some discrepancies at the lowest q values where finite size effects appear to be still relevant. 

IV. CONCLUSIONS 

A noticeable progress in the microscopic study of the ground state structure of medium-heavy, doubly closed shell 
nuclei has been made in the latest years. Going beyond the simple mean field picture and using realistic hamiltonians 
seems now to be within reach of many-body theories. In this context, the variational approach and the correlated basis 
function theory are among the most promising tools. Fermi hypernetted chain and single operator chain techniques 
appear to deal with the complexity of the interaction and of the correlation with the same degree of accuracy that 
has been achieved in nuclear matter. This opens the field both to the microscopic study of dynamical quantitites, as 
the nuclei cross sections, and to the extension of the variationl method to other interesting systems, as, for instance, 
hypernuclei. 
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TABLE I. One- and two-body density sum rules in ^®0 and *''Ca for Jastrow and operatorial correlations. 

Si S2 Sr 

'^O /i 16.00 1.00 -1.00 

fe 16.01 1.05 -0.94 

^Ca /i 40.00 1.00 -1.00 

fe 39.86 1.09 -0.98 



TABLE IL Ground state energies for doubly closed shell nuclei with the S3 potential and Jastrow correlated wave functions. 
The EUL and ACA lines give the energies with r^ dependent and independent correlations, respectively. Energies in MeV. 



(H> E/A (E/A) 



expt 



EUL -33.2 -3.84 

ACA -15.5 -2.36 



«0 EUL -119.7 -8.20 -7.98 

ACA -82.8 -5.89 

^Ca EUL -381.2 -9.78 -8.55 

ACA -262.1 -6.81 



^Ca EUL -394.9 -8.43 -8.67 

ACA -272.3 -5.87 



^Pb EUL -1761.4 -8.50 -7.87 

ACA -1056.3 -5.11 



TABLE III. ^''O and ''"Ca ground state energies per nucleon and radii for the Argonne U14 + Urbana VII and the Argonne 
«8 + Urbana IX models with the FHNC/SOC, Cluster Monte Carlo (CMC) and Coupled Cluster (CC) methods. 

^^O E/A (MeV) rms (fm) 

A14+UVII 

1^0 FHNC -5.97 2.44 

CCM -6.90 2.43 

CC -6.10 2.86 

A8'+UIX 

^'^O FHNC -5.41 2.67 

expt -7.98 2.73 

'*°Ca FHNC -6.64 3.39 

expt -8.55 3.48 



FIG. 1. CBF charge densities and two-body distribution functions for the A8'+UIX interaction. 
FIG. 2. CBF and experimental Coulomb sums for finite nuclei and nuclear matter. 
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